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Michell Trusses Study, Rice U. summer 2013- Summary 

The Michell Truss problem concerns the nature of a structure, constructed 
of bars and cables, of least cost that will balance a given set of point forces. 
The present subcollection of modules resulted from work of two NSF- 
VIGRE groups at Rice University. VIGRE is a program of Vertically 
Integrated Grants for Research and Education in the Mathematical Sciences 
funded by the National Science Foundation. The first was from 2008-9 
academic year research PFUG*, who produced the Connexion module 
Michell Trusses. The second group, consisting of 4 undergraduates and a 
professor, worked during June, July 2013. They studied in more detail 
various aspects of problems related to Michell trusses, made some slight 
corrections to the previous module under the name Derived copy of Michell 
Trusses and wrote five additional modules. *A PFUG is a group of 
Postdocs, Faculty, Undergraduates and Graduate students formed around 
the study of a common problem. 


The connection modules in this collection are: M47207, M47208, M47209, 
M47210, M47223, M47228. 


Derived copy of Michell Trusses 

This report summarizes work done as part of the Michell Trusses PFUG under Rice University's VIGRE 
program. VIGRE is a program of Vertically Integrated Grants for Research and Education in the 
Mathematical Sciences under the direction of the National Science Foundation. A PFUG is a group of 
Postdocs, Faculty, Undergraduates and Graduate students formed around the study of a common problem. 
In this PFUG, we study cost minimizing trusses. We ask what is the structure, constructed of bars and 
cables, of least cost that will withstand a set of point forces. This work was studied in the Rice University 
VIGRE class MATH499 in the Spring and Fall of 2008 and by the 2013 Michell research group at Rice 
University 


Introduction 


The Michell Truss PFUG studies a variational problem posed by mechanical engineer Anthony George 
Maldon Michell in the early part of the last century; what configuration of bars and cables is needed to 
withstand a system of balanced point forces is most economical? Suppose that some collection of forces 
with their points of application (point forces) is given and you are asked to build a structure which can 
withstand these forces. Your choice of materials are bars and cables. Bars can withstand compression 
(forces parallel and pointing inward) but will break under extension (forces parallel and pointing 
outward). Cables can withstand extension but will buckle under compression. You may also choose the 
strength of the bar or cable; the strength is the largest compression or extension the bar or cable 
respectively can withstand. A structure built of such cables and bars will be called a truss. 


The first question the PFUG has addressed is whether one can build such a truss at all. Notice that a bar 
or cable by itself can withstand only one set of point forces; that of two opposite forces applied to the 
endpoints of the bar or cable. These point forces have the property that if the points of application were 
attached to a stationary body, the center of mass would not accelerate and the body would not rotate. This 
property is called balanced. Since point forces can be linearly superimposed, a necessary condition that a 
set of point forces be withstood by a truss is that it be balanced. Below we show that being balanced is 
also a sufficient condition; if the set of point forces is balanced then there is a truss withstanding it. 


Next, the PFUG formulated necessary conditions for a truss to be economical. The cost of a bar or cable 
is its strength times its length. The cost of a truss is defined to be the sum of the cost of bars and cables 
constituted by the truss. In general, there are many trusses which withstand a given a set of point forces. 
The question is, which is most economical, i.e. which truss costs the least? In general this is a very 
difficult question to answer because a minimizing sequence, which we know to exist from the first part, 
may converge in a larger class of measures than described here and it is difficult to formulate any class of 
perturbations of a truss. Below we describe one class of perturbations on trusses with corners which 
yields a surprising necessary condition for economy. A set of perturbations to be studied in the future are 
also described. 


Note:The 2013 Michell group has made some minor corrections and additions to the original module 
created by the 2008 PFUG, but this module is still largely the work of the 2008 PFUG. 


Notation 


We will restrict ourselves to points and vectors in three dimensional Euclidean space, R*. A point force is 
a vector valued measure vd(a) where v is a vector and a is a point in Euclidean space and 6(-) is the 
Kronecker delta mass. If y is a continuous vector field on Euclidean space we define 


Equation: 


(vd(a),) = v(a)-v 


and extend this definition linearly to a sum of point forces. A beam B is a pair of distinct points a and b 
in Euclidean space with a weight w € IR. We can associate with B a mass 
Equation: 


Cost(B) = |w||a — }| 


and point forces field 
Equation: 


a—b 
6B = w——(6(a) — 6(b)). 
cay (Ca) 80) 
6B corresponds to the reaction force of a cable if w < 0 and of a bar if w > 0 with endpoints at a and b. 
Note that (5B, y) = 0 if y is a constant vector field. To see this, note that 


Equation: 
(5B,¢) = (w= 5 (a) -w— 5@),9) 
= 9 (a) oP =F — 9) we 
= 0 


because y(a) = y(b) since » is a constant vector field. We choose the notation 6B to be consistent with 
the notion of first variation of mass, although in this case it differs by the sign of w. A truss T is a finite 
collection of beams and we define Cost(T’) and dT by extending the definitions [link] and [link] linearly; 
Equation: 


Cost (T) = S> Cost (B) 


BeT 


and 
Equation: 


67 = S° OB. 


BeT 


If f is a point force field (a sum of point forces) and T is a truss, then T is said to equilibrate f if 6T = f 
in the sense that 
Equation: 


(67, y) = (£, 9) 


for all continuous vector fields y.f = fd (a1) +--+ + f,,0 (a,) is said to be balanced if 
Equation: 


b 
Vo fi =0 
i=l 


Equation: 


ae xa, = 0. 
i=1 


Existence of equilibrating trusses 


One may easily check that 6B is equilibrated if B is a beam and by linearity 6T is equilibrated if T is a 
truss. The first natural question we deal with is the converse question; is every balanced point force field 
equal to dT for some truss T’? In other words, can any balanced point force field be equilibrated by a 
truss? This section answers this question in the affirmative. The sufficiency will follow from a proof by 
induction. 


Lemma 1 Let f = fd(a) + gd(b) be balanced. Then f is equilibrated by a truss T. 
Let T consist of the single beam (a, b) with w = g- (a — b)/|a — 6. Then [link] implies f = —g and 


then [link] implies that f is parallel to a — b, because f x (a — b) = 0. 
g:(a—b) f-(a—b) _ |-flla—b| 


0S aa ea | f| because f is parallel to a — b. Therefore f = cr because 
Waa is a unit vector parallel to f. Thus f = we By definition 
Equation: 
a—b 
6T = w (5 (a) — 6(b)) = fd (a) — fb (b) = fd (a) + gb (0) =F. 


ja — | 


Clearly then [link] holds. 


Lemma 2 Let f = fd(a) + gd(b) + hd(c) be balanced with a, b and c not lying on the same line. Then f 
is equilibrated by a truss 7’. 


Without loss of generality, a, b, and c lie in the zy-plane and a = 0.f, g, and A must lie in the plane, 
because otherwise it would be impossible to balance the forces. Because a, b, and c do not all lie on the 
same line, two of the points, say 6 and c, are linearly independent; dotting [link] with b and then c implies 
that f -e; = g-e3 = h-e3 = 0 where e; is the unit basis vector parallel with the z-axis. Hence f can be 
expressed as a linear combination of b and c 

Equation: 


b Cc 
f =w— +u4.—. 
|| Ic| 


Consider the point force field e = e,d (b) + e,d (c) where e, = g + we Ty ande, =h+ wat We claim 


e is balanced; 


Equation: 


beg = Gt hwy 


and 
Equation: 


ep Xb+e.Xc = (stove) x b+ (a+ w.& )x 


b 
=gxb+ bx uy thxetex wer 
Cc 


= 9gxb+hxe 
=0 


because 6 is parallel to Wo BT (and likewise for c) and f is balanced. According to Lemma 1, there is a 


truss R with equilibrating e. Let S be the truss consisting of the collection of beams (0, 6) and (0, c) with 
weights w, and w, respectively. We claim T = R+ S equilibrates f. Let y be a continuous vector field. 
Then 


Equation: 
6T =6R+6S 
=e + y-(5(6) — 6(0)) + we (6(c) - 8(0)) 


[b| 


5(c) — 6(0) (5 tue) 


= 95(b) - wrrgp8 (0) +S (6) ~ we 5(0) + on 75 (0) + we 5(0) ~ £6(0) 
= f5(a) + 96(b) + hd(c) 
=f 


which then implies [link]. 
The following lemma shows when the volume of the parallelpiped spanned by three vectors can be made 


to be nonzero by sheering the parallelpiped in a fixed direction by an appropriate amount. An example in 
two dimensions is drawn below. 


at+ty 


b+n: 


Lemma 3 Let a, b, c be distinct, nonzero points in Euclidean space and v a nonzero vector. Then there 
exists ¢ € R such that v € Span{a + tv, b+ tu, c+ tv}. 


If v € Span{a, b, c}, then just take ¢ = 0. If any of the vectors, say a and 8, are linearly dependent, then 
for some constant, n, b = na. For any nonzero t, a + tv is linearly independent from b + tv because 

b+ tv = an + tv is not equal to n(a + tv) = an + ntv) since ntv ¥ tv because v and t are nonzero. 
Applying this argument again to b and c if necessary, we can get three vectors, a + tv,b+ tv, andc+ tv 
, which are guaranteed to be linearly independent. Because the span of any three linearly independent 
vectors is all of R%, every vector is contained in Span{a, b, c}, so it must be that v € Span{a, b, c}, as 
desired. 


Combining the above lemmas, we have 


Theorem 1 A necessary and sufficient condition that a point force field f = fd (a1) +--- + fd (a,) 
be balanced is that it is equilibrated by a truss 7’. 


(Necessity) Suppose f is equilibrated by the truss T. Let vy € R® and y(a) = v forall x € R®. Then 
[link] implies 
Equation: 


Bt Bt 


0=(6T,y)=(f,9)=S fi-vav- So fi. 


i=l i=1 


because y is a constant vector field. Since v was arbitrary, it can be chosen to be 
Equation: 


yw 
v= Dh 


so that 
Equation: 


and conclude that 
Equation: 


b 
ee 
t=1 


which implies [link]. For w € R?, define the skew symmetric matrix 
Equation: 


i W3 0 Wi 
—W2 WI 0 


It is easy to check that tv = w x v. Let yy € R° and y (x) = Zv. Then 


Equation: 
‘ . a; — b; ' ~ Ff a; — b; 
(6T, ) = Y= (6B, ¢) = So wi (y (ai) — g (0:))- ——T = So wi (a b;)y- a 
i=1 i=1 la; — i=1 ja; — B;| 
Le . 
= So ui (ai ip 2G 
= ja; — b5| 


because the cross product w; (a; — b;) x v is perpendicular to the vector a; — b;, which is parallel to 
a;—; 
|a;—bj| 
Equation: 


so their dot product is zero. By [link], this implies 


Le 


lL lL le 
0 = (6T, ¢) = (f, 9) =Vofi-v(a)=S fi -Guv= So fia; xv= Sofi X A,:V. 
i=l i=1 i=1 


i=1 


Since v is arbitrary, we infer [link]. Thus f is balanced. 


(Sufficiency) We proceed by induction on v. If v = 2, then f is equilibrated by the truss 7’ guaranteed in 
Lemma 1. If vy = 3 and ay, a2 and a3 do not all lie on a line, then f is equilibrated by the truss T 
guaranteed in Lemma 2. If vy = 3 and aj, ag and a3 lie on a line, consider the equivalent point force field 
g = fid (a1) + fod (a2) + fd (a3) + 06 (a4) where az is a arbitrary point chosen off the line 
containing a1, a2 and a3. The result for v = 4 proves that g is equilibrated by a truss JT’. However, since 
(g, y) = (£, v) for all continuous vector fields y, this implies that T'’ also equilibrates f. Assume that 
sufficiency holds for v > 3. If f,,1 = 0, then we are done. Otherwise, according to Lemma 3, there is 

p € Rso that 

Equation: 


frsi € Span (a1 — (av41 + pfr41), a2 — (av41 + pfr+1),a3 — (Qr41 + pfr41)). 


Let b = a,41 + pf,4, and let 
Equation: 


fry =); (a, = b) + 25 (a _ b) + 23 (a3 _ b). 


Ay+ + pvr 1 


vel 


Define the point force field g = gi6 (b1) + --- + g 6 (b,) by 
Equation: 


%: = f,+ (a; —)2,, a= 1,2,3, 
9: = fi, t=A,-++ yy, b; = a,, t=1,-+-yv. 


We claim that g is balanced; by [link] 


Equation: 
Vy Vy 3 Vv 
SO Gi = f+ > (Gi — 12 =o fit fu =0 
i=l i=l i=l i=l 
and by [link] 
Equation: 


Vy Vy 3 
So 94 Xb; = So fix ait S52; (a; — b) x a; 
i=1 i=1 i=1 
Vy 3 3 
= So fix ai + > Qj (a; — b) x (ai — b) + $2; (a; — b) xb 
i=1 i=1 i=1 


= Sofi x ay + Fea x (p44 + pfr41) = 0. 
i=1 


By induction, there exists a truss R equilibrating g. Let S be the truss consisting of the collection of 
beams a, — b, ay — b, a3 — b and b — a,,,, with weights 2,, 22, 23 and |f,,| resp. Arguing as in Lemma 
2, we find that T = R + S equilibrates f. 


Some Economical Trusses 


After proving the question of existence we turn to the question of economy. We say that a truss T is 
economical if it satisfies 
Equation: 


Cost(T’) < Cost(S) 


whenever 6S' = dT. That is to say, the cost of T is less than or equal to that of any truss which 
equilibrates the same force system as T’. We begin by describing some global statements about 
economical trusses that can proven by choosing special test vector fields y in [link]. Then we make local 
perturbations on trusses with corners to find a necessary condition for economy. 


Some Global Properties 
A surprising and easily proven fact is 


Lemma 4 A truss is economical if the weights of the beams are all of the same sign. Such a truss lies in 
the closed convex hull of the support of the point forces it equilibrates. 


To prove the first statement, note that 
Equation: 


(57, x) =) (6B,2) = yw (ad) = Dvla~ 


BeT BeT BeT 
= S> Cost (B) = Cost (T). 
BeT 


if w > 0 for each B € T. On the other hand, if S is a truss with 6S’ = dT, then 
Equation: 


Cost (T) = (6T, x) = (65, x) = Sw |a — b| < Cost ($). 
Bes 


Hence T is economical. 


To prove the latter statement, let K be the closure of the convex hull of the support of the point forces 
equilibrated by 7’. Then Kk is a convex polyhedron; let H be the hyper-plane passing through one of its 
sides. Without loss of generality, assume that H is the zy-plane and 

H, = {x € R®: 2; = x2 =0,23 > 0} is the upper half space, and Support (6T) C R* \ H,. Let 
yp (x) = e3 ® e3x for 3 > 0 and 0 otherwise where eg is the unit basis vector (0,0, 1) and e3 @ eg is 
the tensor product 


Equation: 
000 
e3@e3= 0 0 0 
001 
Then, 


Equation: 


0 = (52, y) = Sow eg @e3(a—B) = Yow 


BeT 


which implies that w = 0 if it corresponds to a beam lying in H,. 


Some Local Properties: Cutting Corners 


Our object of interest are two dimensional trusses with corners. We have shown that an economical truss 
cannot have comers. To show this, we analyzed a perturbation of a truss with corners which consists of 
cutting the corner and replacing it with a flat top. Specifically, we showed that for any corner, the cut can 
be made sufficiently small so that the perturbed truss costs less. This surprising result suggests that any 
economical truss, if made of both cables and bars, has as boundary a differentiable curve and is supported 
on a Set of positive two dimensional area. 


We define a comer to be the union of three beams {B,, Bz, B3} C T which share an endpoint p, lie ina 
halfplane about p and (dB; + d6Bz + 6B3, yp) = 0 for all y with Support(y) Cc U for some 
neighborhood U of p. 


By rescaling, translating and rotating, we may assume that B, and B2 form two sides of an isosceles 
triangle and one endpoint of Bs lies in the base of this triangle. The base of the triangle and B, form an 
angle of 8 and the base of the triangle and B3 form an angle of y. The height of the triangle is /. We will 
need the relation 

Equation: 


Fer O>0. 


The sum of B,, Bz and B3 will be called the four point truss, 74. The four point truss equilibrates the 
system of forces f1, f2 and f3 with points of application the intersection of B,, Bz and B3 resp. with the 
base of the triangle. We require 

Equation: 


Ji hijo tis = 0 


or 
Equation: 


\fi |+| fo |=|fs |sin yg esc 8, | fi |—|fo |=|f3|cos yp sec 0. 


The simplicity of this truss allows an for a straight forward calculation of the cost. For 2 = 1, 2, 3 let w; 
and A; denote the strength and length resp. of beam B;. It is easy to see that for this truss to equilibrate 
the the strengths of each of individual beam must be equal to the force that is based on a point of the 
beam and aligned to it, thus: 

Equation: 


w1=|fil, we=|fel, ws = —|fsl- 


The lengths of the the Ty are also easily calculated. B; and Bz being the equal sides of an isosceles 
triangle it might be superfluous to state that: 
Equation: 


Ai = A2 
With the use of trigonometry it is clear that 
Equation: 
Ai = Aq = I\csc O| 
and 
Equation: 
A3 = I\csc »| 


The cost of Ty is given by the formula 
Equation: 


3 
Cost(T1) = S> |wilA; 
i=1 
=|f1 |l csc 6+] fe |l csc O+|f3|1 csc y 


=l\f3 (sin y csc? 6+ esc y). 


We have used [link] in relating | f;| and | f2| to | fz|. Note that the cost is linear in J and diverges as 9 + 0 
ory — 0. 


Four Point Truss Five Point Truss 


Having calculated the cost of the T4 we perturb the corner in order to produce a structure that may yield a 
lower cost compared to that of the 74. Our main objective being that of investigation of the consequences 
or eliminating corners from an already existing truss. 


The modified truss will be obtained by removing the corner of the T4. A horizontal line parallel to the 
triangle base line a distance / — h from the base of the isosceles traingle is added. The beams B‘, and Bi, 
are formed by shortening B, and Bz to where they intersect this added line. To mantain a truss structure 
B3 must be replaced by three new beams. B, connects B’ and Bz in the place where the cut of the 
corner was made. B% shares a connections point with B, and Bi while B3 is connected to the node that 
is shared with both B4 and B). The shape of the new truss that is obtained is that of a trapezoid, the 
applied forces remain same as for T4 We shall hence forth refer to this truss structure and the five point 
tent truss, Tsp. 


The complicated nature of the 75, makes the calculation of the lengths 4,'s and strengths w,'s slightly 
less obvious than they were for the original 74. However, easy trigonometry and some calculations yield 
both the length's and the strengths of each individual beam. We use the same parameters to describe this 
truss as before, namely y, 8 and J. However, we need three more parameters to calculate the cost. These 
parameters are the following: h which describes the height of the cut and y and 2 represent the angle 
between the base and By and B3 resp. 


Keeping in mind some of the properties that govern the truss we have calculated the length's of each 
individual beams to be: 
Equation: 

A = (l—h) csc 6 

AL = (l—h) csc 8 


Xs, = y/ (Leot g — h cot 6)? + (1h)? 


Mt = yf (cot o + hot 6)? + (0h)? 
A4 = 2h cot 0 


The magnitude of strengths are similarly calculated to be: 
Equation: 


wil= | fi 

wo |= | fal 
ws|= | f,| sin 8 csc py 
w3|= |fo| sin 0 csc p_ 


|wa|= |fi| cos 0+ |fi| sin 6 cot py, 


Having acquired all the information we require we can now calculate the cost of the T5 ;: 
Equation: 


Vv 


Cost (Ts,n) = So wir; 
i=1 
= |fi| (l—h) csc 6+ |fo| (J — h) csc 0 
+|f1| sin 8 csc eV cot y — h cot 6)” + (1— hy’ 


+|f2| sin 6 csc p_- cot y + h cot 6)? + (l- hy? 
+|f1| cos 6+ |fi| sin 6 cot yp, (2h cot 6) 


In order to check that the perturbation lowers the cost for some sufficiently small h, we differentiate this 
expression and evaluate the derivate at h = 0. Even further simplification reduces the formula to the 
following somewhat complicated expression: 

Equation: 


d 
qh cost (Ts,n) n-o = |fil[— csc 6+ sin 6 cos” ~ (cot 6— cot ~) 


— sin O(cot y cot 6+ 1) +2 cot A(cos 6— sin 6 cot »)| 
+|fo|[— esc 6— sin 0 cos” y (cot y+ cot 6) 
+ sin (cot y cot 6 — 1)| 


With the help of the above analsysis, we have shown that 
Equation: 


d 
an oo (T5,n) n-o0 < 0 


independent of J, and y. This proves the claim that two dimensional trusses cannot have corners. 


Other Deformations 


Other deformations of trusses worth considering are perturbations of other junctions. Consider the case 
when B,,---,B,, are beams meeting at a point p. Suppose we shorten each beam by introducing a 
polyhedron about p whose m vertices meet the beams B,,---,B,, and each edge represents the location 
of a new beam. From Lemma 4 above we see that the cost does not change with the perturbation if the 
junction consists of only cables of only bars. This study would help rule out certain cases when both 
cables and bars meet at a junction. 


A second deformation consists of making piecewise affine deformations within a triangle as follows; 
divide the triangle into three sub-triangles with common vertex in the center. Choose a base and let A, be 
the affine transformation of R? which maps points x to points A p (£) whose distance to the line passing 
through the base is scaled by p. The affine transformations in the remaining triangles is uniquely 
determined by continuity. Thus, 7, (a) be the unique piecewise affine transformation of R? with 

Np (a) = a for z not in the triangle and 7, («) = A, (a) for z in the first sub-triangle. Applying 7, to a 
truss and then adding the necessary beams to ensure that is equilibrated then yields the perturbation. 


Tensegrity 


Tensegrity and Trusses 


Introduction 


The term “tensegrity”, coined by Buckminster Fuller in the 1960s, is a combination of the words 
“tension” and “integrity”. In a tensegrity system, the compression members do not touch each other, and 
are stabilized by the tension members. Skelton and Oliveira, in their paper [link] define a class-k 
tensegrity system as “a tensegrity system such that no more than k rigid bodies make connections at a 
given node (with frictionless ball joints).” 


Due to the high strength and low weight of tensegrity structures, Fuller imagined building 
“unprecedentedly large” structures [link], perhaps even large enough to cover a city. Kenneth Snelson, a 
former student of Fuller's, is renowned for his tensegrity sculptures, which are on display around the 
world. In addition to manufactured structures, tensegrity systems are abundant in nature. The elbow is a 
class 3 tensegrity joint, while the foot and shoulder are class 2 tensegrity joints. In a red blood cell, the 
protofilament is the compression member and the spectrin dimer are the tensile member [link]. 


This module further examines properties of tensegrity systems by looking at the paper “Optimal 
tensegrity structures in bending: The discrete Michell truss” by Robert E. Skelton and Mauricio C. de 
Oliveira. The following analysis and diagrams come largely from their paper, with changes made by the 
2013 Michell group to help clarify their ideas for readers at the undergraduate level. 


Michell Spirals 


Consider a set of line segments, connected as shown below. In the diagram, y = 7g and 6 = =. 


Definition 1 Let r, define a set of radii from a common origin, 0, for 2 = 0,1, 2,...,q. Let pe, 
£=0,1,2,..,q—1 define the lengths of lines beginning at points with radius r¢ and terminating at 
points with radius rz,,;. Then a Michell spiral of order q is defined by the connections of lines of length 
pe, satisfying 

Equation: 


Tei =are, Pe=CPre, £=1,2)2049; 


where a,c > 0. 


By a visual inspection, one can observe that 
Equation: 


Te41 COS p+ py cos B= Te 
and 
Equation: 


rpiz Sin yp = py sin B 


where y is the angle between each radius and £ is the angle between the beam and the radius. 


Therefore, using some algebra and trigonometry, one can get that 
Equation: 


sin B sin y 
a= ———— and c= 


sin (6 + ¢) sin (6 + y) 


In order for the spirals to converge to the origin, it must be that a < 1, which corresponds to 
y+28 <n. 


Michell Topology 


The magnitude of a node is determined by the radius on which it lies. That is, 
Equation: 


|| mae | =|] amen || 


forali+k=m-+n. 


Notice from the diagram above that nodes with the same radius are related by a phase-shift of 2my, 
where m is an integer. Using complex notation, this can be expressed as 
Equation: 


j2 
Ni+m,k—m = e? OP nay 
The line segments comprising the spiral are described by defining the vector connecting nodes n;z and 


Nik+1), 
Equation: 


Mik = Nik — Nik+1 
where the vector n;, has magnitude and phase given by 
Equation: 


Nik = niger, Nik =Titk, Prix = (a ~~ k)p 
such that rz satisfies the conditions in the definition for some specified rg. The mirror image of all lines 
reflected about the horizontal axis are obtained by computing the conjugate of the vectors m,,, and 


shown in [link] below. 


Note that, by a visual inspection, mj, = pipe? I6+(—k)¢] | where Pm, m= 0,1,2,...,q, satisfy [link]. 


Definition 2 A Michell Topology of order g consists of Michell spirals of order less than or equal to q 
and their conjugate spirals where [link] and [link] hold. 


Force equilibrium at a generic node 


A force applied to a node can be represented by the vector w;, applied to the node n,, and written in 
complex form as 
Equation: 


Wik = Wiperrvs 
where 
Equation: 


Pwix = Dik + (i — k)p 


and 0;x is the angle at which the external force wx is applied, measured from the radial line to node nx. 


At node njx, the forces along directions mj;,, My;, mM; 4-1, Mz,;-1, respectively, have magnitudes 
firs tri, fi,r—1, thi—1. Therefore, the sum of forces at node nj; yields 
Equation: 


=j2-17 — 


— tei-1 T 


M ik 


Multiplying [link] by the vector p;,,e4(*-9” gives 
Equation: 


Ditk fine” +tget? = fps —tyg eh) + we | = 0. 


Because both the real and imaginary parts of the above equation must be zero, 
Equation: 


| cos B cos | @ a | cos(8+y) cos(6+ | Lo tek [vs mae 7 (3) 
—sinf sin B! \ fix BaP | git (6+) sin(6+y)} \fin-1 Pitk | in 6x! “ Oo 
Linear propogation of forces 


Define a vector x, € R2(*+)) which contains forces in all members that lie within the radii r, and Tot 
. That is, 


toe 
Sao 
ta-1 


I cect 


—— 
bia—i 


fous 


tao 


foe 
Define a vector Uy such that 


Wa+1,0 
Wa,1 
Wa-1,2 
Uy = Wa—2,3 


Wa-3,4 


W0,a+1 


Theorem 1 Let a truss be arranged according to the Michell topology of order q, satisfying the 
conditions in the definitions, having external forces w;, applied at the nodes n;, with 7, k > 0 and 
i+k < q. Let x, contain the forces(normalized by the member length) in all members within the band 
of members between radii r, and r,,,. Then the forces propagate from one band to the next according 
to the linear recursive equation Xgi1 = AgXq + Bug for some matrices Ag and By. 


Michell topologies under a single bending load 


If tixfik < 0, then one member of the truss is in compression while the other is in tension, so the truss is 
said to be bending. 


Theorem 2 Let a truss be arranged according to the Michell topology of order q satisfying [link] and 
with a and c given by [link]. Assume the only force applied on the truss is w = we”? at node ngp. If 
Equation: 


p>0, 0<B<|6l<x-B, yp+2B<7, 
then the forces at all nodes satisfy 
Equation: 
tinfix <0 forall i<q, k<i 
Furthermore, if 
Equation: 
B<0<n-8 
then t;, > 0 for allz < q, k < 7. Otherwise, if 
Equation: 


B-17<0<-8 


then t;, < Oforalli<q,k <i. 


First, take the case where gq = 1. In this case, noo is the only node at which external forces are applied: 


In the case when 0 < 6 < 4, there are three cases of interest. The first is when |Ao0|< 6B < 4. Then 
sin (28) > 0 and 
Equation: 


—1t << —268< O09 —7n<0 ==> sin (000 — 8) < 0 
Equation: 


0<O0+8<28B<07r => sin (O00 + 8) > 0 


which implies tog < 0 and foo < 0, or, in other words, that both members of the truss are in tension. 


The second case is when |Ao0 — | < $< 4.A similar analysis shows that to9 > 0 and foo > 0 so 
both members of the truss are in compression. 


Lastly, when 8 <|6o9|< a — 8 then too foo < 0, indicating that one member of the truss is in tension 
and the other in compression, so the truss is said to be bending. 


Now, consider the case where gq > 1 and assume that 0 < 8 < 099 < a — f'so that tog > O and fog < 0 
(the opposite case is handled similarly). Consider the node ng;. In the force diagram, wo, = 0 by 
assumption, hence ¢1,-1 = 0 because ng is at the lower boundary of the truss. Repeat the analysis from 
before for the fictitious force Wo = | foo|moo. Because the truss has a < 1,8 < 09, = p+ 6 <7, s0 


that t19 > O and fo; < 0. The same idea, applied to all nodes on the upper boundary with k > q, proves 
that tx9 > O and fox < 0 for all k < g. 


The analysis at node nyo is similar. Repeat the above analysis for a fictitious force Wi9 = — |too|Moo. 
The vector has a negative sign because the angle 049 is measured clockwise, so the vector needs to be 

flipped to get the same angle @ as before. Therefore, for all nodes njg with  < q, tox > O and fin < 0 
forall k <q. 


Continuing on, for the node n4j, use a fictitious force wi; = —|to1|Mo1 + | f10|mm10. Because a < 1, 
one can conclude that 8 < 6;; < a — to get that t3; > O and fy; < 0, and continuing in this fashion 
gives t;; > Oand fi; < 0 forallz < q. 


The last type of node to consider is an interior node of the form nj; where i # k # 0. At the node njz 
create the fictitious force wiz = —|ti-1,x|/M20 + | fii{maii. The assumption that a < 1 leads to the 
conclusion that 8 < Oj, < a — 8, yielding t,; > O and fj, < 0. Therefore, when 8 < 099 < a — B, 
Equation: 


te >O and fix <0 forall i,k<q 


A similar procedure for —@ > 009 > 8 — 7 yields 
Equation: 


te <0 and fx >O forall iik<q 


Material volume of optimal structures 


Let o;, 7 =1,...,m, denote the yield stress of the material used to construct the 7th member, m,;, 
loaded with force density 7; and with constant cross section area A ;. To avoid material failure, the force 


density at any given member should satisfy 
Equation: 


with equality holding in minimal volume structures. Some algebra yields 
Equation: 
vj =|| m; || Aj = |? 


1 . 
——o;(|| my ||) |] mj ||", J=1,..m 
OF 


where v; is the volume of each member. Therefore, the total volume of a structure, V, is 
Equation: 


Let b; and s; be the bar and string vectors for the 7th compressive or tensile member in a structure with 
my bars and m, strings. Therefore V = V; + V;, where 


Equation: 
Mp 1 ' Ms 1 5 
Vs = $5 =Aj(|l Bj I) |] by I? and Ve = $2 (| 5 II) I ss | 
j=l Aj j=l V5 


Assuming that all bars are the same material, dj = ) forall j = 1,..., mp. Likewise, if all strings are 
made of the same material, then Y; = 7¥ for all 7 = 1, ..., ms. Hence 


Equation: 
™mob 1 Ms 1 5 
Vo = D5 =Aj(Il By II) by I? and Vs = $7 Syl 55 Il) Il s5 Il 
j=l A j=l v 
so that 
Equation: 
Mp 1 ; Ms 1 5 
V=Ve+Ve= >) AGI By II) Il By I? +] alll ss II) Il's5 Il 
jal Aj jal VS 


Now, let J be defined by 
Equation: 


J =2\7V + (7 = ) y n: wi 
=1 


From [link] and [link], it is true for a structure in static equilibrium that 
Equation: 


n Ms mp 
doa wi = Yo vill $5 IN) Mh sal? — So AC By IN) | Bs I? 
=1 j=l j=l 


Let B= $75", Aj(I| by |I) || b; ||? and S = 30%" -y;(Il 8; ||) || sz ||? Therefore 
Equation: 
_—1 1 = 
J = 27] B+ =s| + (7-2) [Ss - Bl 
nN Y 

= 27B4+2\8 +79 —7B-AS+ XB 

=7B+AS+7S5 +.B 

= (X+7) (B+S) 


Thus 
Equation: 


J= (X+7) 


do valll $4.11) Mss? +95 As bs ID) I by | 
j=l j=l 


Assuming that \ = 7, we get that 
Equation: 


f= IV 


™Mp 1 Ms 1 
do Aull Bs I) I Bs 1? +2 Sas(ll 83 II) Il ss | 
=e, j=l A 


= 2 bs Aj(I] By Il) I By II? + 5 ra(ll 85 I) Il sy | 
j=l j=l 


= 2)? 


which is consistent with [link]. This implies that minimization of J is independent of the choice of 
material properties, and depends only on the variable topology. 


The Steiner Tree Problem 
An overview of the Steiner tree problem. 


Introduction 


Given a set of points, called terminals, a spanning tree is the set of terminals 
plus a number of edges, in which two terminals connected by exactly one 
simple path. The problem of the minimal spanning tree is to find the shortest 
spanning tree, where the length of a tree is the sum of the lengths of its edges. 
In some cases, points in addition to the original set, called Steiner points, can 
be carefully placed so that the spanning tree connecting the new Steiner 
points and the original set of points is shorter than any tree that could be 
created by adding, removing, or moving any Steiner point. The resulting 
spanning tree is called a Steiner tree. 


The Michell truss problem begins with a set of point forces, which can be 
viewed as points in the Steiner tree problem. Because shorter beams cost less 
than longer beams of the same weight, it is beneficial to minimize the length 
of the beams used in a truss. Applications of the Steiner tree problem may 
help optimize a truss by minimizing the length of its beams. 


Definition 


The Toricelli point of a triangle, also called the Fermat point, is found by 
constructing two equilateral triangles on any three sides of the triangle, and 
drawing a line from the new vertex of each equilateral triangle to the opposite 
vertex of the original triangle. The location where these two lines intersect is 
the Toricelli point. 


Each triangle has exactly one Toricelli point. In an acute triangle, the Toricelli 
point is inside the triangle, and the lines from the Toricelli point to each 
vertex of the triangle create 120° angles. Toricelli proved the shortest path 
connecting all three vertices of a triangle consists of three lines, each from the 
Toricelli point of the triangle to a vertex of the triangle. 


Theorem 1 


No two edges of a Steiner tree can meet at an angle less than 120°. 


Proof 

If the two edges meet at an angle of less than 120°, then select points on the 
edges that are equidistant from the vertex, and locate a new Steiner point at 
the Torricelli point of the triangle created by the vertex and the two selected 
points. Because the triangle is acute, the three new edges created by 


connecting the Toricelli point to the vertices of the triangle must be shorter in 
sum than the total length of the original two edges. 


Theorem 2 


A Steiner tree has no crossing edges. 


Proof 


If a Steiner tree had crossing edges, the cross would create at least two angles 
less than or equal to 90°. 


Definition 


The degree of a point in a tree is defined as the number of edges connected to 
that point. 


Theorem 3 


Each Steiner point of a Steiner tree is of degree 3. 


Proof 


From Theorem 1, no two edges meet at an angle of less than 120°, so the 
Steiner point can be of degree 3 at most. A Steiner point of degree 1 can 
simply be eliminated, as it creates an unnecessary edge, and a Steiner point of 
degree 2 can just be replaced with a straight line. Therefore a Steiner point 
must be of degree 3. 


Theorem 4 


A Steiner tree connecting n terminals contains at most n — 2 Steiner points. 


Proof 


The Euler characteristic for lines is 1, which means that 

#vertices — #edges = 1. Therefore there are n + k — 1 edges in a tree, 
where n is the number of terminal points and k is the number of Steiner 
points. 


Each Steiner point has 3 edges coming from it, and each terminal point has at 


least one, so there are at least Bhd total edges; the division by 2 accounts for 
the fact that each edge is counted at two vertices. Therefore 
Equation: 


Fete me pee 


which, by simple algebra, yields 
Equation: 


n—-2>k. 


Definition 


A Steiner tree with n — 2 Steiner points is a full Steiner tree. 


Theorem 5 


Every Steiner tree that is not a full Steiner tree can be decomposed into a 
union of full Steiner trees. 


Proof 


See [link] 


Examples 


In a 3 by 4 rectangle, two Steiner trees can be created, but the one with edge 
between the Steiner points parallel to the long side is the Steiner minimal tree. 


In a 1 by 2 rectangle, the minimum spanning tree is not a Steiner tree. 


Definition 


A topology on a graph determines how many terminals and Steiner points 
there are, and how the points are connected to each other. A topology is a 
degeneracy of another if it can be obtained from shrinking edges of the 
original. The set of degeneracies of a topology G is denoted by D(G). 


A topology is a Steiner topology (full Steiner topology) if every Steiner point 
has degree 3 (and terminal has degree 1). 


If F is a full Steiner topology, then D, (F’) is the set of Steiner topologies in 
D(F). 


It is easy to see that a Steiner tree contains no closed shapes, called cycles. 
Otherwise, one edge could be deleted, and the resulting structure would be of 
shorter total length while still connecting all of the points. 


Theorem 5 


A Steiner tree whose topology is in D(F’) for some full Steiner topology F is 
the unique minimum tree among trees whose topologies are in D, (F). 


Proof 


Let T be a tree whose topology G € D, (F). Let v1,..., vu, be an 
enumeration of the vertices. Define f;; = 1 if [w:, V5 is an edge in G and 
fi; = 0 otherwise. Then the length of T’ is 

Equation: 


f= Shale y 


i<j 


|v; — v,| is a norm, and norms are convex due to the triangle inequality, so T’ 
is strictly convex except when edges remain the same direction ina 
perturbation. In that case, when the three edges from a Steiner point are 
parallel, moving the Steiner point towards the side of the two overlapping 
edges decreases the total length, and thus |T'| is strictly convex. Because a 
Steiner tree is a local minimum, strict convexity guarantees that it is the 
unique minimum among trees with topologies in D, (F). 


Corollary 


There exists at most one relatively minimal tree for a given topology. 


Definition 


A Steiner hull for a given set of points N is a region which is known to 
contain a Steiner minimal tree. 


In the Euclidean plane, the Steiner hull is bounded. 


Lemma 


Let |v;, v;| be any path in a Steiner minimal tree, where v; and v; are defined 
as before. Let L(v;,v;) be the region consisting of all points p satisfying 
Equation: 


|p—v; |<|v; — vj| and |p — vj |<|u; — 9; 


L(v;, v;) is the lune-shaped intersection of circles of radius |v; — v,| centered 
on vu; and v;. No other vertex of the Steiner minimal tree can lie in L(v;, v;). 


Proof 


If q were such a vertex, the Steiner minimal tree would contain either a edge 
from q to v; not containing v;, or vice versa. If the SMT contains a edge from 
q to v;, for example, the SMT can be shortened by deleting |v;, v] and adding 
[g, vj], a contradiction. 


Similar statements can be made to further specify where a Steiner hull may 
be. 


Remark 


It is easy to see that a full Steiner topology with n terminals has 2n — 3 
edges. Let f(n), > 2, denote the number of full Steiner topologies with 

mn — 2 Steiner points. Adding a Steiner point to the middle of any of the edges 
and connecting it to a new terminal results in a full Steiner tree for n + 1 
terminals. Geometrically, this causes the edge to become two separate edges, 
connected by the new Steiner point. In order to be consistent with the fact that 
edges must meet at 120° angles at Steiner points, the two new edges will 
rotate to meet at an angle of 120° with each other and with the edge 
connecting the new Steiner point to the new terminal. Thus there are 2n — 3 
ways to create a new full Steiner tree, so 

Equation: 


f(n + 1) = (2n — 3) f(n) 


which has the solution 
Equation: 


Let F'(n, k&) denote the number of Steiner topologies with n terminals and k 
Steiner points such that no terminal is of degree 3. Then F'(n, &) can be 
obtained from f(k) by first selecting & + 2 terminals and a full Steiner 
topology on it, and then adding the remaining n — k — 2 terminals one at a 
time at interior points of some edges, creating a "kink" in each edge. The first 
terminal can go to one of the 2k + 1 edges, the second to one of 2k + 2 


edges, and so on, until the (n — k — 2)th point is added to one of the 
2k+(n—k—2)=n+k-— 2 edges. Thus 


Equation: 
Pe ae (n+k— 2)! 


Letting n3 denote the number of terminals of degree 3 
Equation: 


F(n) aS (2 )F (n= nak + ng) SEE 


A table containing the values of f(n) and F(n) for n = 2, ..., 8 is given. 


n 2/3 | 4 5 6 7 8 
f(n) 1]),1]3 15 105 945 10395 


F(n) 1 4 31 360 5625 110880 2643795 


Computational Complexity 

The optimization problem: 

Given: A set NV of terminals in the Euclidean plane 
Find: A Steiner tree of shortest length spanning N 


can be recast as a decision problem: 


Given: A set NV of terminals in the Euclidean plane and an integer B 
Decide: Is there a Steiner tree T' that spans N such that |T'] < B 


The Steiner tree problem is NP-complete. This means that it is at least as hard 
as any problem in NP, but the solution is still verifiable in polynomial time. 
However, the problem cannot be solved in polynomial time. The discrete 
Euclidean Steiner problem, however, is not known to be NP-complete. 


Given: A set NV of terminals with integer coordinates in the Euclidean plane 
and an integer B. 


Decide: Is there a Steiner tree T' that spans NV such that all Steiner points have 
integer coordinates and the discrete length of 7’ is less than or equal to B, 
where the discrete length of each edge of T’ is the smallest integer not less 
than the length of that edge. 


There is a version of this problem that is in P: Given a finite set in a Banach- 
Minkowski plane, a minimal Steiner tree for this set can be found in 
polynomially bounded time with an algorithm that executes only graph 
theoretic functions. 


The Steiner Ratio 


The Steiner ratio is the smallest possible ratio between the total length of a 
minimum Steiner tree and the total length of a minimum spanning tree. It is 


conjectured that, in the Euclidean Steiner problem, the Steiner ratio is * 


This is called the Gilbert-Pollak conjecture, and was though to be proven in 
1990, but some people do not accept this proof. Gilbert and Pollak also 
conjectured the ratio for higher dimensional spaces, but this was disproven. 


It's easy to see that this is true in the four point case. Alternatively, the Steiner 
ratio is sometimes given as the total length of the minimum spanning tree 
over the length of the minimum Steiner tree. The Steiner ratio for the 
rectilinear Steiner tree, which is when all edges of the Steiner tree are 

3 


perpendicular, is >. 


Other applications 


Steiner trees have many biological applications. They can be used to describe 
the ways proteins connect and fold. Steiner trees are also used in phylogeny. 
Organic chemistry diagrams often include Steiner trees. For example, CHgN3 
is a regular hexagon Steiner tree. 


Steiner trees have many computer science applications as well, such as in 
network design and circuit layout. 


A Look at a Paper by Strang and Kohn 

This module attempts to present the paper "Hencky-Prandtl Nets and Constrained Michell 
Trusses" by Gilbert Strang and Robert Kohn in a way that is accessible to undergraduate 
mathematics students. 


Hencky-Prandtl Nets and Constrained Michell Trusses 


Introduction 


The paper “Hencky-Prandtl Nets and Constrained Michell Trusses” by Gilbert Strang and 
Robert Kohn [link] provides many interesting insights into the problem posed by Michell in his 
1904 paper “The limits of economy of material in framed-structures” [link]. This module 
attempts to present the former paper in a manner that can be understood by undergraduate 
students in mathematics. 


Michell's problem assumes a material with given tensile and compressive yield limits, say a9, 
and allows concentrated forces. A force F' can only be carried by beams with cross-sectional area 
F'/oo, so if the points of application of the forces are too close together, there will be difficulty 
in finding enough space. Because there is no restriction on the magnitude or the spacing of the 
external loads, which are surface tractions f distributed along the boundary I’ of the given 
simply connected domain £2, one needs to be careful to ensure Michell's problem is self- 
consistent. To do this, we remove the limitation imposed by 09 and allow increasingly strong 
beams whose cost is proportional to their strength. This way, the cross-sectional area doesn't get 
too large, so there isn't a problem with finding enough space. 


Variational forms of the stress problems 


The following is a brief derivation of the Michell problem, where the stress constrains come 
from the equilibrium equations: 


Equation: 
divo=0in Qio:n=f on TI 
where 
Equation: 
o= 
Tay Oy 
Equation: 


and its eigenvalues, or principal stresses, are 


Equation: 


1 1/2 
A12 = . (0, Gyo ((c2 — Oy) + 47” zy) ) 


The formula for [link] comes from the quadratic formula for the 2-D Lagrangian strain tensor. 
Bars are placed in the directions of principle stress, which, by definition, come from the 
orthogonal eigenvectors of o. The required cross-sectional areas at any point are proportional 
(scaled by a9) to |A;| and |A| and the total volume of the truss is proportional to 

Equation: 


@ (co) = is |+]A2 |) da dy 
2 


Michell's optimal design can be determined from the solution to 
Equation: 


Minimize (oc) 


subject to divo=0,0:n=f 


To express the problem in a way that suggests a computational algorithm, we introduce a stress 
function (2, y). For any statically admissible stress tensor o there is a function w such that 


Equation: 
oe ( Vy a.) 
—Way Wax 


where the subscripts indicate partial derivatives. It is easy to see that div o = 0. In the first 
column, for example, 
Equation: 

0 O 0 _ oO 
T dy 


( Woy) = Dyya — Voss = 0. 


w is determined up to a linear function, because all second derivatives of a linear function are 
zero. Given the unit normal n = (nz, n,), the unit tangent t = (—n,, nz), and the condition 


YyyNa — VryNy = fi or (grad w,) -t = fi and similarly for fz. Therefore the tangential 
derivative of wy, is, by definition, Py =f ory P) = f a fi ds by the fundamental theorem 
of calculus, and similarly for f2. In order for the constraints to be compatible, the external loads 
fi and f2 must be self-equilibrating, and therefore the boundary values 7, and 7, are well 
defined: the integrals around the closed curve I’ are zero. Provided that there is no net torque 


from the external forces, ~ itself can be determined. Write ~ = g, WW, = h, where w,, is the 
normal derivative of 2. 


If the signs of A, and 2 are the same, then 
Equation: 


[Ar |+]A2 |=|A1 + A2 |=|Pex + Pyy| 


because the trace of the matrix o is equal to the sum of the eigenvalues. 


If the A; and Ag are of opposite signs, then 
Equation: 


Aa] + |A2| = |A1 — Ag! 


1 
= E(t =s (oz —o,)? + 4r2,) 


1 z 
= 5|¢ ((c. _ Oy)? + 47>, ) 
i 
a ((c. a Gy) as 4r2,) ; 


Therefore, in order to minimize the largest principle stress, the Michell problem becomes 
Equation: 


Minimize / max(|Pz2 + Py 
2 
subject to Y= G9,¥n =honl 


1/2 
, (Wer — Py)? + 402,) ) da: dy 


In real life, unlike in Michell's problem, beams in a truss cannot be arbitrarily strong. Thus, if A 
denotes the upper bound on the stresses, the design problem for a constrained Michell truss is 
Equation: 


Minimize / / rx (2) |+]A2 (2)| 
‘e 


subject to |Ai|< A, divo=0, o:n=f on IT 


The stress tensor is made up of two tensors: the mean hydrostatic stress tensor and the stress 
deviator tensor. The hydrostatic stress is the average of the normal stress components, oz and oy. 


The deviatoric stress is what's left over, namely c? = o — $ (a, + o,)I. The constraint on the 


strength of the beams can act only on a” instead of o, therefore ignoring the hydrostatic 
component of the force. This is useful for working with a composite material that can only 
withstand pure pressure. The constraint then applies to the eigenvalues of o?. Using the same 
quadratic formula as in [link], we get that 

Equation: 


2 


A (0)| =|r2 (0)? = (Fe. “)) rea 


The matrix 0” has trace zero, so the eigenvalues are of equal magnitude and opposite sign. 
Therefore the plain strain redesign problem is 
Equation: 


Minimize // Mt (o”) |+|A2 (o”) 
2 
subject to divo=O0in 2,0:n=f onT, |r (0°)| = )A2 (o”) |" < k? 


Rewriting this in terms of the stress function w, and denoting it by || w ||, we get 
Equation: 


Iv? = (Fle —Ym)) +95 


Therefore the constrained problem above takes the form 
Equation: 


Minimize 2 || || dx dy 
if 


subject to || ¢v ||< kin Q, P=g, Vn=honl. 


If the k is too small, or f; and f are too large, then the material may be too weak to withstand 
the load, and no truss can exist without collapsing. 


The dual problems and optimality conditions 


Now we give an informal derivation of the dual problem, which is a maximization over 
displacements u = (1 (x, y), u2 (x, y)). No boundary conditions are imposed, because 
ao : n= f was prescribed on all of I”. To solve for the maximum over the constraint div o = 0, 


we introduce a Lagrange multiplier. This is a system of two equations, we iF uae = 0 and 
ray a ae = 0, and we denote the multipliers by u, and w2. Then 
Equation: 
min // Ai |+|A2| =minmax // Ay |+|A2|+u- div o 
div o=0 on=f wu 
o:n=f 


=maxmin as |+|A2|— (e (w),o)+ fu-f ds 


F : . Ou; F 
where ¢ is the deformation tensor with components €;; = > ( On + 3) and the inner product 


(¢,0) is }) >) €:;04; = trace (eo). 


The next step is minimization over a, so € remains fixed. 
Equation: 


min |), (0)|+|Xo (2)| — (€,0) 


First, we must diagonalize ¢ by a principle axis transformation, ¢ = R~‘e'R where R-! = R¢, 
so that the eigenvalues €1 and €2 (the principal strains) appear in the diagonal matrix e’. The 
inner product (€, 0) and the eigenvalues \(o) are invariant if this transformation is applied at the 
same time to a: 

Equation: 


(€,0) = trace (eo) = trace (R-te’/RR™'o'R) = trace (e’0’) = (e', 0’) 
Thus the minimization simplifies to 
Equation: 
min|A; (o’) |+|A2 (0')|-e104, — £209 


If je1|> 1, take 0}, = Ke, and 04, = 0 for some K > 0 so that the expression now becomes 


Kle1|—K]e1 is < 0, and let kK — oo so that the minimum goes to —oo. Similarly, the minimum 
is —oo if |e2|> 1. Otherwise, if |e1 |, |e2|< 1, the minimum is zero. To see this, let o;; = 0 if 
le] < Lando}, = Ke; if |e;|= 1. 


Substituting in this minimum, the dual problem is now 


Equation: 


Maximize / (uifi + uefe) ds subject to |e; (u)| <1 
r 


It's expected that this maximum value is equal to the minimum of [link], corresponding to the 
minimum volume of the truss. 


Consider the case when |e;|= 1. When €; = 1 and €2 = —1 the principal strains have equal 
magnitude and opposite sign, which can happen only for a special class of displacement fields. 
There is a condition sot = 0 on the angle 6 between the horizontal and the direction of the 


: Ou Ou : 
strain €; = 1, and the secondary property €; + €2 = 3? 4 By = 0. Here, the geometrical 
problem is identical to that of slip lines in plane strain, where the eigenvalues of stress equal 
tog. In the cases when €; = €2 = 1 or €y = €2 = —1 there is a similar correspondence with 


pure hydrostatic pressure in the stress case. 


Because [link] is equal to 0, for every admissible o and u, 


Equation: 
min // dif Aal=max fu f ds 
div o=0 u r 
o:n= 2 

so that 

Equation: 


[fis (o)|+\A2(0)| de dy > fw- fas 


with equality holding if the optimality condition is satisfied. 


If o and e(w) are simultaneously diagonal, the possibilities for o can be restated as 
Equation: 


if A, (oc) =0, then |e; (u)| <1 
if A;(0) #0, then e;(u) = sgn d; 
because in a diagonal matrix, A; (0) = 0%;,;. 


The dual of the problem [link] becomes 
Equation: 


A1 (2) |+|A2 (2)| — (€, 0) 


min 
|A;|<A 


and the principal axis transformation leads to 
Equation: 
. / | 
yee At (’) |4/A2 (0') |-e1041 — €20%9 


Because 4, cannot be arbitrarily large, —oo is no longer reached when |e;|> 1. Hence, there are 
now three cases, depending on the value of |e;|: 
Equation: 


if |e;| >1, then A; = A sgn e; 
if |e;| =1, then A; =|A,| sgn €;,0 < |A;| < A 
if je;| <1, then A; =0 


In the second and third cases, the minimum value in [link] and [link] is still zero. However, in the 
first case, when |e;|> 1, the minimum is no longer —oo. When o’ is diagonal with 

o,, = A; = A sgn «;, |A;|—E;0); = A(1—|e; |). Thus, the minimum in [link] and [link] can be 
expressed as 

Equation: 


Ma = Amin (0, 1—|e, |) + A min (0, 1—|e2 |) 


Here, A; is being chosen through the minimum expression depending on the fixed ¢;. Hence the 
dual problem can now be stated as 
Equation: 


Maximize [fms dx ay+ | (uifi + uo fe) ds 
yi 
\0) 


The cases of interest arise when one family of bars is in tension and the other compression, 
meaning A;A2 < 0. Then there is a combination of Hencky-Prandtl nets, one coming from slip 
lines, and the other from Michell trusses. 


The slip line net will occur when €; > 1, €2 < —1, A1 = A, Aq = —A. In this case, || w || 
represents the difference s (Az — Az), which has the constant value A, which is the yield surface 
in plane strain. The Michell situation occurs when €; = 1, €2 = —1, A1 = 0, Az = 0. 


The other constrained problem [link] differs from the dual only in the condition of 
incompressibility: 


Equation: 


Maximize [fm dx ay + f (ui fi + U2 fe) ds 
6 A 
subject to div u=0 


The kinematic constraint div wu = 0 arises because the trace of o is unconstrained in [link]. The 
optimality conditions relate the principal strains €; to the eigenvalues Ae =X; ie ) :e anda? 
are simultaneously diagonal, and 


Equation: 
if le| > 1; then ba =k sgn €;, 
if |e| <1, then A? =0. 

Open problems 


e Find the most economical truss given a balanced point force field. 
¢ Does the stress field in an optimal structure vary smoothly, so that the lines of principal 
stress are C’! curves? The problem of regularity runs deeper than that of existence. 


Michell Truss Stress Measure 

More detailed and accessible derivations of the equations in section 2 of the paper 
"Michelle Trusses and Lines of Principal Action" are provided here. See Stress and 
Measure Theory modules for background material. 


Michelle Trusses 


Unidirectional Curve Stress 


Let C be a simple curve in R®. In other words C is the image of a map 
r € C1 ([0, 1], R2) such that # (t) A 0 for each t € (0, 1] and r is injective. We will 
also require without loss of generality that r(t) trace out the curve C at a constant 
1,. 
speed, making |7 (¢)| constant on [0, 1]. Therefore |7 (t)| = Junie = #1 (C) for 


each t € [0,1] and the unit tangent vector to C at time t will be 
Equation: 


The curvature « at time t is then given by 
Equation: 


si dr/dt _ #(t)/HE1(C) _ FL) | 
|dr /dt| KH" (C) (#1 (C))? 


We now define a measure o that is proportional to the stress along the curve C’. 
Equation: 


Thus o© is a rank 1 symmetric matrix measure. 


We will now treat o@ as the functional on C, (R?, R**4) that corresponds to the 
measure o@. For any € € C, (R?, Res). 
Equation: 
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To find the divergence of the measure 0 


ueCi (R?, R**4) be given and we evaluate the functional —div o 
Equation: 


, we let an arbitrary vector-valued function 
“ [ul 


(—div a) [u) = o° [Vu] 


Vudo"l 
Re 


| 


es rae (Wu (r (t));# (t) @ #(t)) at 


Using the definitions of the inner product and the tensor product, we can convert the 
matrix product (Vu (r (t)); 7 (t) © 7 (€)) into a vector product: 
Equation: 


(Vu (r (#))# (4) @F(E)) = (Vuln (t))s# ()F()” ) 


Therefore 
Equation: 


(-aiv o°) uy) = —~— [ (Ser )r@)) at 
WHC) I) \at 


= [burda fi tury aia— fo (usm) ase ie 
= 7dplul — 7d, [ul] —K.7} lc [ul 


(r5B — 764 —-K.%0} lc) [24], 


where A = r(0) and B = r(1). Since this equation holds for any 
ueC} (R4 ; Rese), we say that the measure —div o° is given by 
Equation: 


—div o° = 76g — 764 —K.%' | ¢. 
When the curve C is just equal to the line segment [A, B], 7 = Beal and « = 0. 


Therefore our equation becomes 
Equation: 


[link] 


Stress in a Truss 


The stress in a beam is given by +a/4-#! for some . A finite collection of such 
8 Y > 


beams is a truss. Therefore the stress in a truss that consists of | nodes at 
{Aj, Ag, ; a , A;} is given by 
Equation: 


where we require that Amm = 0 for each m. 


Recall that the stress o of a structure that is in equilibrium must satisfy the restrictions 
—div o7 = Fando =o’. Itis easy to see that o will be symmetric by definition. 
Therefore, if our truss is in equilibrium, there must be a force distribution F’ such that 
Equation: 


l 
F = —div o = —div S° = Ann glAm An] 


m,n=1 


Since the operator —div is linear, we get 
Equation: 


l 
fe = S- Amn div ol Amn] 


m,n=1 
l 
Amn —A 
= 5 (nr pe eae 
dy, 4m Ay) |Am — An] 


which will be our new requirement for a truss to be balanced [Link]. 


Cost of a Truss 


If we are given a force distribution F’, a set of points 2& = {A,, Ao,---,A;}, anda 
set of weights A = {Amn : m,n = 1,---1} such that 


P= ae Amn (64,, — 94, ) ane then the cost of the truss structure will be 


the "total mass" of the structure, which is given by fy, |o]. 
Equation: 


lo|(#") 


= “i( U Ar An] 


l<m<nx<l 


= > [ol ((Am; Anl) 
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~ Ds Apg 74041] ([Ams An) 
1<m<n<l|p,q=1 
= Prmnaltet + ram ol4etel| (Am; Anl) 
l<m<n<l 
A, —A A, —A 
— 2? ike Ecie |i aot a os [AgAd ([Am, An]) 
1l<m<n<l |Am = A,| |Am -_ A, 
l 
(An — A,) ® (Am — An)| 1 
= 7 Primal Am Aa) Am = And geet pas] (Any An) 
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m,n=1 
Cost Minimization 
Given a truss T, which consists of a set of points 2a = {A,, Ao,---, Az} anda set of 
weights A = {Amn : m,n = 1,---1}, and given a force distribution F of finite 


support, we define a set pars 2 ) of all admissible stress measures for the truss T. 


Equation: 


Up (2) = {2 CM (2, 5x?) : -div 0 = F, and o = 3 —Amn saute 


m,n=1 


l 


m,n=1 


AmAn 


The requirement o = > Ain o! is equivalent to the restriction that 7 


be rank 1 and supported on a finite collection of simple curves. 


We also define the set 


Equation: 


Dp (2) = {eo & M (2, oo :—div a= Fi, 


of all stress measures for structures that balance the force distribution F’. The stresses 
in this set do not necessarily correspond to trusses, since we have dropped the 
requirement that the stress be rank 1 and concentrated on curves in R@. In fact, the 


stress measures in dp (2) may be spread out over a region of R?. Also, we may 


allow the given force F to be diffuse. 


Our original optimization problem was to find a truss T° such that the cost of T° is 
equal to inf { inf { flcj:c€ OE (@ \ \ for a given force distribution F of finite 


support. This is equivalent to finding a stress measure 0° € LJ (=f (2)) such that 


[|o"\=int i 


Such a 0°, however, does not always exist [link]. 


Equation: 


oO 


reer (a))} 


T 


Tensors, Lagrange methods, and Hausdorff Measure 
Presentation Notes Harry 
The Stress Tensor: 


Informally, tensors are geometric objects that describe linear relations 
between vectors, scalars, and other tensors (say, linear maps or dot 
products). A tensor can be represented by a multi-dimensional array of 
values and the number of indices needed to specify a component in this 
array is said to be the order of the tensor. 


Multi-dimensional arrays: 


1. Scalar [a] - this is a tensor of order 0 because we don't need any indices 
to specify "a" within a 1x1 array 


2. Vector [ai ] i=1,...,n - This is a tensor of order 1 because we need 1 index 
to specify a component of the vector 


3. mxn matrix - this is a tensor of order 2 
4. Tensors of higher order can be thought of as multidimensional "boxes" 


Notice that we can create a tensor of a new order through the multiplication 
of tensors. For example, the products of a nxmxp tensor with an px1 vector 
will yield an nxm matrix. 


Definition: The stress force on an object is said to be the force per unit area 


So, the stress tensor is an order 2 tensor (a matrix) which takes a unit length 
direction vector normal to the surface of an object as an input and gives the 
stress force vector (with respect to the given coordinate system) at this 
surface as an output. 


Then, for our problem, which is in 3 dimensions, the stress tensor is a 3x3 
matrix with entries Tij, where i,j=1,2,3. 


Imagine a beam and say that we would like to calculate the stress vector at a 
point, p, in the beam due to a force at the endpoint. Then, given a 
coordinate system, consider an "epsilon-cube" around p. Given an applied 
force at the endpoint of the beam (which translates to the same force at p) 
Tij is the stress vector acting on the ith face in the -jth direction - note that 
opposite faces have the same index. For an explanation of these properties, 
see [1]. 


One can check that in the context of a Michell Truss, the stress tensor will 
be a rank-one matrix and that the system of forces is in equilibrium if the 
divergence of the stress tensor is the negative of the applied forces. In other 
words, admissible trusses satisfy 

Equation: 


Vo -—-F 


where F is the vector sum of the applied forces. See [5] for an explanation 
and derivation of this equation. 


Definition: Strain is the amount of deformation an object experiences as a 
result of external forces when compared with its original size and shape. 


The strain tensor is operates similarly to the stress tensor, but the entries Tij 
are replaced by partial derivates of the position of a point, p, with respect to 
movement in the coordinate directions as a result of external forces. In the 
context of a Michell Truss, the strain tensor will always be a rank-one 
matrix as well. 


One can investigate the relationship between stress and strain via a 
stress/strain curve. 


Euler-Lagrange Equations 


How do we optimize functionals? Consider a functional 
Equation: 


over a region containing a and b. 


Let u(x) be a minimizer for I(u). Then u(x) should satisfy the Euler- 
Lagrange equations as follows: 


Equation: 
d df(zx) ! GINS) is u’ (x 
ds | an’ Gey ROM OD) | = Gacey BO) 909) 


For derivations and explanations of these equations, see [2] and [3]. 


Note that this doesn't mean that any function like u(x) satisfying the 
equations is a minimizer. In fact, any stationary function for I(u) will satisfy 
the Euler-Lagrange equations. However, the point for us is that if we can 
create some sort of functional describing possible truss structures given a 
set of balanced forces, then we can employ the Euler-Lagrange equations to 
see which of these truss structures are "stationary," and thus possibly cost 
minimizers for supporting the given set of point forces. 


The next question one might ask is how to construct such a functional for a 
finite set of point forces. 


Using Lagrange Multipliers to Model a Finite System of Point Forces 


The Lagrange multiplier method is used to find the maximum and minimum 
of a function subject to a set of constraints. For example, one could find the 
maximum and minimum values of the function f(x,y) = 6x + 8y subject to 
the constraint g(x,y) = x*+ y?- 1 = 0 (e.g. maximize and minimize the 
function on the unit circle). Our goal is to minimize the cost of an 
admissible truss structure subject to constraints regarding the application of 
a finite number of point forces. For example, given a balanced system of 
point forces in the plane, 


Unknowns: 

1. Weights of each possible beam in the structure with wij = wji 
2. Li multipliers to balance the x component of the applied forces 
3. Mi multipliers to balance the y component of the applied forces 


Constraints = balance of force components in each coordinate direction at 
each point of application. These equations are of the form: 


At a point, a0: 


(x-component force, y-component force) = Li [sum of the x-components of 
the weights of each beam connected at a0 | + Mi [sum of the y-components 
of the weights of each beam connected at a0 | 


As one might imagine, this problem can become complicated to solve in 
cases which lack symmetry. In the usual basic equilateral triangle example 
(with a point also at the centroid) that we used, there were still 4 constraint 
equations (one for each point of application) with an x-component part, y- 
component part, and all the unknown weights. 


One should note that any solution for the simultaneous equations will 
generate an admissible truss structure. However, solving these equations 
simultaneously is a difficult problem in higher dimensions, or even when 
there are just more than four points and forces, none of which are 
symmetrically arranged. Overall, this is usually an impractical approach for 
finding a cost-minimizing truss structure given a set of balanced point 
forces. 


Intro to Hausdorff Measure 


Definition: If U is a nonempty subset of $%”, the diameter of U is said to be 
|U = sup{ |x — Y|, %&Y € U} 


Definition: Fix 6 >0. If EC UU; where i = 1,...,00 and |U|< 6, we say that 
the collection of U;'s is a 6-cover of E. 


Definition: Let EC R”. For every 6 >0 anda >0, 6,a€ R, let H3(E) = 
inf{S> |U;|%, EC Uj, |U; |< 6} where i = 1,...,00. Then, we get the a- 
dimensional Hausdorff measure of E by letting 6 go to 0 as we take the sum 
over all d-covers. 


So, H°(E) = lims_.9H5‘(E) = sup{ Ho} among all 6 >0 since HY 
increases as 6 decreases. 


One important basic result in developing the notion of Hausdorff measure is 
that if 6 is strictly less than the distance between two positively separated 
sets, E and F, then no set in a 6-cover of EF U F can intersect both E and F. 
This means that for some a >0, 


H¢ (EU F) = He (E) + Hg (F) 


Now, suppose we have two values for a, s and t, with s<t and we want to 
consider the Hausdorff measure of a set EC R”. Then, H(E) >ds — tH} 


(E) = if H3(E) >0, then H£(E) is infinite as 6 goes to 0. 


From this, we can see that there is a unique value for a, called the 
Hausdorff dimension, d, of E such that H°(E) = 00 if 0< a < dimE and 
A°(E) = 0 if dimE < a < o For all of the above in the Hausdorff measure 
section, see [4]. 


So given a set E, how do we calculate its Hausdorff dimension? 


Definition: Let f:A—B (and for our purposes, assume, A,BC R”). The 
If(z)-F)| 


function, f, is said to Holder continuous of exponent ¥ if: sup ray 


Fs 


oo among all x,yeE A. 


Proposition: If f: A—-B is onto and Holder continuous for some ¥, then 
dim(B)< = dim(A), where these are the Hausdorff dimensions of A and B. 


Proof: Let UU;, i = 1,...,00 be a 6-covering of A for some 6 >0. Then, f(U 
U;) covers B and by the Holder continuity condition, we have that for each 
U;, C|AN Uz|">|f(AN U;)|, where the bars denote the diameters of those 
sets and C is a constant. 


This implies that for any a, C7 JAN Uz|°>|f(AN U;)|7 


Taking this sum over all 1, we get that 57| f(A U;)|7<C7 SJANU,|*< 
C7 S\|U;|“>H7(B) <C7 H(A) 


Now, suppose the a is larger than dim(A). This means that H(A) = 0 so 
that H°(B) = 0 as well. Decrease @ until it reaches dim(A). Then, dim(B) 
x dima, as desired. 


We can use these measures in the context of the Michell problem. 
Developing the problem in new ways might allow one to use different 
techniques to find minimal cost structures among admissible configurations. 
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